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Abstract 

We consider the one-dimensional Kardar-Parisi-Zhang (KPZ) equation with half 
Brownian motion initial condition, studied previously through the weakly asymmet- 
ric simple exclusion process. We employ the replica Bethe ansatz and show that 
the generating function of the exponential moments of the height is expressed as 
a Fredholm determinant. From this the height distribution and its asymptotics are 
studied. Furthermore using the replica method we also discuss the multi-point height 
distribution. We find that some nice properties of the deformed Airy functions play 
an important role in the analysis. 

1 Introduction 

Surface growth phenomena appear widely in nature and have attracted much attention 
in non-equilibrium physics. In 1986, Kardar, Parisi and Zhang proposed a stochastic 
differential equation which describes surface growth with local interaction [1]. For one- 
dimensional case, the KPZ equation is given by 

dh{x,t) X fdh{x,t)y d'^h{x,t) 



dt 2 



(^) .^^.v/I^,(M). (1.1) 



Here h{x,t) represents the height of the surface at position s G M and time t > 0. The 
first term represents the effect of nonlinearity and the second one describes the smoothing 
effect of the surface. ri{x,t) represents randomness described by the Gaussian white noise 
with covariance, 

{7]{x, t)r]{x', t')) = 5{x - x')5{t - t'). (1.2) 
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The parameters A, u, D determine their strengths. For the one-dimensional KPZ equa- 
tion (1.1), it has been shown by the dynamical renormalization group that the height 
fluctuation scales as 0{t^^^) as t goes to inflnity [1]. The exponent 1/3 coincides with the 
ones found in Monte Carlo simulations of many stochastic models of surface growth. The 
KPZ equation is accepted as a prototypical equation describing the universality class called 
the KPZ universality class. 

Our understanding of the KPZ universality class has deepened in the last decade. Not 
merely the exponent but the height distribution functions have been computed based on 
intriguing connections with the random matrix theory [2] . For the totally asymmetric sim- 
ple exclusion process(TASEP), which is an exactly solvable model in the KPZ universahty 
class, the current distribution function for the step initial condition has been obtained 
exactly and it was found that in the long time limit, it converges to the GUE (Gaussian 
Unitary Ensemble) Tracy- Widom distribution function [5], the largest eigenvalue distribu- 
tion in random hermitian matrix [2]. 

An interesting finding from the studies of distribution functions is that they detect 
the difference of initial conditions, which the scaling exponent of the height cannot. For 
example, in the surface growth model called polynuclcar growth(PNG) model, it has been 
recognized that the GUE Tracy- Widom distribution describes the height distribution in 
the droplet growth [3] while in the flat initial condition, it is described by the GOE (Gaus- 
sian Orthogonal Ensemble) Tracy- Widom distribution [4] . Furthermore, generahzations of 
these results to multi-point distribution function have also been studied and the universal 
processes, the Airy2 [6] and Airyi [7, 8] processes have been obtained. For recent progresses 
of this topic, see [9, 10, 11]. 

Since 2010, the studies on the KPZ equation and KPZ universality class have entered 
a new stage [12]. First as an experimental progress, high-accuracy measurements for the 
exponents and the height distribution in the KPZ growth problem using turbulent liquid 
crystal have been performed and the GUE Tracy- Widom distribution function was observed 
as well as the critical exponents [13]. Second, we have begun to understand exactly height 
distribution function of the KPZ equation itself. In [14, 15, 16, 17] the height distribution 
at a single position was computed for the narrow- wedge initial condition, 

A/,(^,i^O) = -M, S^O, (1.3) 

from which the surface grows to a parabolic shape. In the long-time limit, the fluctuation 
around the macroscopic shape is described by the GUE Tracy- Widom distribution. The 
analysis of [14, 15, 16, 17] is based on recent progress on the current distribution of the 
(partially) ASEP [18, 19] and a fact that, in the weak asymmetry limit, the stochastic time 
evolution of the current of the ASEP can be mapped to that of the height described by 
the KPZ equation [20]. For more recent developments, see [21, 22, 23, 24, 25, 26]. 

More direct approach without relying on the results on the ASEP has also been devel- 
oped recently [27, 28, 29]. The idea was first proposed in [30], in which the author used 
the interesting relation that the exponential moment of the height is represented as the 
dynamics of the one-dimensional (^-function Bose gas with attractive interaction, which 
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is solved by the Bethe ansatz [31, 32]. In the original work [30], only the ground state 
contribution was considered and the dynamical exponent 1/3 was obtained. In [28, 29], 
the authors succeeded in taking into account the whole contribution of the eigenstates and 
obtained the Fredholm determinant representation of the generating function of the expo- 
nential moment for the narrow wedge initial condition (1.3). More recently a multi-point 
generating function is also discussed [33, 34]. 

The replica method is quite attractive since we expect to be able to calculate various 
quantities for the one-dimensional KPZ equation directly and easily. In particular, we 
expect it would be a powerful tool to understand the universality in the KPZ equation. 
There are two recent progresses in this direction: In [35], the authors discussed a pro- 
cess characterizing the renormalization fixed point of the KPZ universality class and the 
Fredholm determinant expression of its transition probability has been obtained. In [36], 
the authors have obtained the exact height distribution for the fiat initial condition and 
clarified its convergence to the GOE Tracy- Widom distribution in the long time limit. 

In this paper, we pursue the potential of the replica method by applying it to the half 
Brownian- motion initial condition depicted in Fig. 1. This is one of the typical spatially 
extended initial conditions and is written in terms of the single valued function h{x, t) as 



Here B[x) represents the one- dimensional standard Brownian motion with B{0) = and 
a — (2i/)~^/^AD^/^. For this initial condition, the macroscopic shape expected by solv- 
ing (1.1) without the noise term is 



The main interest in this paper is the fluctuation around this macroscopic shape. We 
expect that a one- dimensional Brownian motion describes the fluctuation in the positive 
region {x > 0) since it is known to be a stationary measure of the KPZ equation. On 
the other hand, in the negative region {x < 0) where the parabolic growth is observed the 
situation is the same as the narrow wedge case. Thus the fluctuation property around the 
origin shows a crossover behavior between the two typical growths. 

This initial condition has already been considered in [38]. There the analysis is based 
on the result for the ASEP with step Bernoulli initial condition [37] and naturally leads 
to an expression for the height distribution in the form of contour integral. We take 
the replica method and mainly treat the generating function of the exponential moment. 
Accordingly our formula is somewhat different from the one in [38] and is expressed in 
terms only of real quantities. Furthermore the replica approach allows us to discuss the 
multi-point distribution function with the help of the factorization assumption by [33, 34] 
and properties of the deformed Airy functions discussed in Appendix. 

This paper is arranged as follows. In the next section, we state our main results. Briefly 
explaining the relation between the KPZ equation and the (5-Bose gas in Sec. 3, we give 




(1.4) 




(1.5) 
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h{x,0) 



Figure 1: Half Brownian motion initial condition 



a derivation of our main result, the Fredholm determinant expression of the generating 
function (Theorem 1 in Sec. 2) in Sec. 4. In Sec. 5, we discuss another expression of 
Theorem 1 stated as Proposition 2 in Sec. 2, which is useful for the compact expression of 
the height distribution function described by Theorem 3 in Sec. 2. In Sec. 6, we discuss 
the multi-point height distribution. In Sec. 7, we consider the interpretation of our result 
as free energy distribution of a directed polymer in random media. Concluding remarks 
are given in the last section. 



2 Model and main results 
2.1 The generating function 

The KPZ equation (1.1) is in fact ill-defined as it is. As time goes on the height profile 
approaches the stationary one described by the Brownian motion, for which the non-linear 
term in the KPZ equation is not well-defined. A proper prescription is proposed in [20], 
which we follow here. In this scheme, one defines the height of the KPZ equation by 

K,\,d{x, t) = {Z^^x,d{x, t)) , (2.1) 

using the solution Zy^\^r){x,t) of the stochastic heat equation, 

Ql = + ^Vix, t)Z^,xM^, t), (2.2) 

which is a well-defined stochastic partial differential equation of Ito-type. This is called 
the Cole-Hopf solution of the KPZ equation because the equation (2.2) is related to the 
KPZ equation through the (inverse of) the Cole-Hopf transformation (2.1). We will give 
more explanations in subsection 3.1 below. Hereafter we investigate the properties of these 
regularized quantities ht,^x^c){x,t) and Z^^x^D{x,t) with the initial conditions (1.4). In the 
following we denote by (■ ■ ■) the average over both ri{x,t) in the KPZ equation (1.1) and 
B{x) in the initial condition (1.4). 
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We are interested in the distribution of the height h,^,x,D{x,t). It is well established, 
known as the KPZ scaling, that the fluctuation of the height scales like 0{t^^^) and non- 
trivial correlations are seen in the x direction with scale 0(t^/^) when t is large. Let us 
define a parameter 7^ which scales as 0{t^^^) and a rescaled space coordinate X by 

7t = {a'ut)K (2.3) 
with a given below (1.4). We introduce the scaled height ht{X) by 

^K,x,D = -^^-ItX^ + lMn (2.4) 

Here the second term corresponds to the macroscopic shape in (1.5); it is an interesting 
aspect of the KPZ equation that one has also to take into account the first term to focus 
on the height fluctuations. 

To study the distribution, it is often useful to consider the generating function of the 
moments {Zi,^x,d{x, t)^), = 0, 1, 2, • • • . We define G^^{s; X) as 

^) = E^=^ (^"^« f?llij)\e-i----' = (e-.-'— (2.5) 
Ar=o ■ \ ^ ^ ' 

By using the replica method, one can express the moment {Zj,^x^r>{x,t)^) of our problem 
in the language of the 5-function Bose gas, which is a well-known exactly solvable model. 
One can further perform the summation over N to obtain the Fredholm determinant rep- 
resentation of the generating function. 

Theorem 1 G^^{s]X) is expressed as the Fredholm determinant with the kernel acting on 



G^,{s-X) = det (1 - PoKxPo) ■ (2.6) 
Here Ps represents the projection onto (s, 00) and the kernel of Kx is given by 



Kx{ujj,UJk) = y(-l)"-' / ^g-n(a.,+a.,)-2iq(a;,-a;,)-7fng2+#n3-7*n. V 2! ^ , 

tl J^-^on ^ r (ig - f + f ) 

(2.7) 

where T{x) is the gamma function and Cn satisfies Cn > X/^t + n/2. 

Here for an operator on L^(R) with kernel K{x,y), the Fredholm determinant is defined 

by 

det {1-K) = J2 ^TTT / ■ ■ ■ / det {K{x„Xk))^,^^ (2.8) 

M=0 



Kx{ij,ik)= / dyAf (^. + y,-,--]Air{ik + y,-,--] _^.„ , ■ (2.10) 



where the right hand side is assumed to converge. 

A derivation of this result will be given in Sec. 4 after the explanation of the relation of 
iVth moment of Z^^x^D^x^t) with particle dynamics of the (5-function Bose gas in Sec. 3. 
For the narrow wedge case (1.3) where Zj,^A,D(a^? 0) = 5{x), the corresponding generating 
function has been obtained in the same form as (2.6) [28, 29]. But the kernel Kx{iOj,uJk) 
is replaced by 

Xnw(^^j,u;fe) = ^(-1)"-' / ^e-"(^^•+'^*)-'^«(^^-^*)-^*'^'''+l^"'-^*"^ (2.9) 

n=l ^ 

which is independent of the position X. 

One can replace the kernel by the one written in the form of products of two deformed 
Airy functions. 

Proposition 2 In (2.6), the kernel PqKxPq can be replaced by PqKxPq where 

1 X\ , / 1 A:\ e^*2/ 
— , Air Ufc + I/, — , -^Z^r— -y^s 

Here Ai^ (a, b,c), Air{a,b,c) are defined by 

Af{a,b,c)^— [ dze'^'+'^riibz + c), (2.11) 

1 C°° 1 

Air(a, b,c)^- dze-+-^ (2.12) 
27r i_<^ r {-tbz + c) 

/n (2.11), represents the contour from — oo to oo and, along the way, passing below the 
pole Zp = ic/b. 

The functions Ai'"(a,6, c), Air(a, 6, c) have appeared in [38]. Note that these become 
the ordinary Airy function if the gamma function factors in the integrand are eliminated 
in (2.11) and (2.12). For the narrow wedge case, the corresponding kernel is the one where 
Ai^(a;, l/'jt, —X/^t) and Air(a;, l/7t, —X/^t) are replaced by the Airy function. To get the 
kernel in proposition 2 from that in Theorem 1, one has to find some generalizations of 
formulas utilized in [28, 29]. They and a derivation of (2.10) will be given in Sec. 5. Some 
properties of the deformed Airy functions are summarized in Appendix. 

2.2 The height distribution function 

The information of all the moments is enough to extract that of the probability distribu- 
tion function. Applying the discussions in [29, 33], we can find a formula of the height 
distribution from the generating function G^^(s,X). Let F^^{s;X) be 

F^,(s;X) = Prob (^^Kx,t) + y| + 7t^' ^ = P^ob(/it < s). (2.13) 
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By using the Predholm determinant (2.6), F^^{s;X) can be expressed as follows [29, 33]. 

/oo 
due-^'''''-'''g,,{u-X). (2.14) 
-oo 

Here 

g^X^^- X)^^. (det(l - PoK+Po) - det(l - PoK^Po)) (2.15) 

where K^(x,y) is the kernel (2.7) or (2.10) in which the term is replaced by — ± ie 
with e > being infinitesimal. 

Using (2.10), and the relation l/{x ± ie) = V{l/x) =F m5{x), where V denotes the 
Cauchy principal value, we can easily find K"^ in (2.15) is represented as 

I X\ .. 1 X\ f 1 



/oo 
dykf ( + -, -- ) Air ( a + -, -- 
-oo 



It It) \ '7*' It) vi-ei'*(" y'^ 



.r / . 1 X\ ( ^ 1 X 



T iTrAi^ + -, Air Ufe + • (2-16) 

V It It) \ It It) 

Substituting this expression to (2.15) and using basic properties of determinant, we even- 
tually arrive at the expression in terms of the Fredholm determinant. 

Theorem 3 

/oo 
due-^'"''-'''g^,{u;X). (2.17) 
-oo 

Here g^^{u]X) is expressed as a difference of two Fredholm determinants, 

g^^{u; X) = det (l - P^iB^ - PX.)P„) - det (l - P„<P„) , (2.18) 

where 



.(^1,6)= r 

^0 



BliCi, 6) = / dyAf (Ci + y, -, --) Air f 6 + y, -, -- 

^ ' ^ It It) \ It It 



Jo ei*y - 1 V V It It) V It It. 

Ai^ fei - --^ Air ^6 - y, -- n , (2.19) 

(2.20) 



It It J \ It It 



PUCi,^2) = Ai^ U, --) Air U2, -, 

V It It) V 7t It) 



In [38], another expression using a contour integral was obtained whereas ours (2.17) arc 
represented as the convolution with the Gumbel distribution. For numerical analysis. 



7 




Figure 2: The probability density functions at the origin X = 0. The three cases for finite 
^'s (7t = 1,2,5) are depicted as the dashed hues while the solid line corresponds to the 
7t oo limit (see (2.22)). 

this form would be convenient. Actually, from this expression, we can readily draw the 
picture of the probability density function dF^^{s] X) / ds as in Fig. 2. In this figure, the 
graphs are drawn approximating the Fredholm determinant by a finite dimensional matrix 
determinant using a simple discretization. For more precise estimation, the method in [39, 
40] is available and actually it was applied to the narrow wedge case [41]. 

The distribution function (2.17) has a similar form as the narrow wedge case obtained 
in [14, 15, 16, 17]: If we replace the functions Ai'"(a;) and Air(a;) by the ordinary Airy 
function in (2.19) and (2.20), the distribution function (2.17) becomes the one for the 
narrow wedge initial condition. 

We also consider the long-time limit {t — )■ oo) of the distribution function. Let us 
remember a basic fact that the probability distribution function is in general written as 
the expectation, Prob(X < s) = E(9(s — X)), where Q{y) is the step function, Q{y) : = 
1 (y > 0), {y < 0). Noticing lim exp[-e"°''] = 9(x), we take the t oo limit in (2.5) 

a—^oo 

and see 

lim ProhCht < s) = lim G^,{s;X). (2.21) 

7t->oo 7t— i-oo 

Combining this relation with Theorem 1 and Proposition 2, one obtains 

lim PToh{ht <s) = det (1 - PslCxPs) , (2.22) 

7t^oo 

the kernel is given by 

^x{ij,ik) ■■= lim Kxi^j - s,^k- s) 

7i— )-oo 

POO 

= / dyAii^, + y)Aii^k + y) 
Jo 

+ Ai(e,) (^e-^+^«^ - 2" dye-^'y Ai{^, + y)^ . (2.23) 
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This distribution function was obtained in the study on the PNG model with an external 
source [42, 43] and the TASEP with step Bernoulli initial condition [44]. At X = 0, another 
expression using the solution to the Painleve differential equation was given in [45]. In Fig. 
2, we illustrate the picture of the probability density function at the origin X = as a 
solid curve. The derivation of (2.22) will be given in Sec. 6 including the discussion of the 
multi-point height distribution. 



2.3 Multi-point distribution function 

The replica approach allows us to discuss multi-point height distribution function. We 
introduce the n-point generating function 

G..(Mn,W„) = (e-^?-^^"'"^^'-^^'), (2.24) 



where we abbreviated Si, • • • , s„ and Xi, - ■ ■ , X^ as and {X}n respectively and we 
set < ^2 < • • • < Xn- 

In [33, 34], the authors proposed the "factorization assumption" for the n-point gen- 
erating function for the narrow wedge initial condition. In this paper, we show that if we 
employ the same approximation to the case of half Brownian motion initial condition, we 
obtain the following result: 

Gl{{s}n,{X}„)^det{l-Q), (2.25) 
and the kernel Q{x, y) is given by 

/oo 
du2--- dur,{ui\e^^^-^^^"\u2) ■ ■ ■ K|e(^"-^i)^Li|ii„+i)$({ii - s}„), 
-oo 

(2.26) 

where H is the Airy Hamiltonian H — — ^ + u, and 

Lj{x,y)^ f dwAir(w + x,-,-^]Af(w + y,-,-^), (2.28) 
Jo \ It It J \ It It J 

with j > 1. In (2.25), wc put the symbol jj in order to represent explicitly the fact that this 
is an expression after using the factorization approximation. The factor {x\e^^^~^''^^\y) 
in (2.26) can be expressed as 

(x]e(^^-^^)^|y) = r dze-'' Air ( x + z, Af ( y + z, , (2.29) 
J-oo V It It ) \ It It ) 
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see (A.l). For the one-point case n—1, however, the method in [33, 34] is exact. Actually, 
when n = 1, wc find the Frcdholm determinant is equivalent to the kernel (2.10): In the 
one-point case, the kernel (2.26) becomes 

Q{ui, U2) = Li{ui, U2)^ui - si). (2.30) 

We divide it into Q = Q1Q2, where 

Q,{u, w) = $(m - s)Ai^ (w + u, ) xoiw), (2.31) 

V It It J 

Q2{w, u) = Xo(u')Air (w + u, -, -— ) , (2.32) 

V It It J 

where Xs{C) — ^ {s < ^ < 00), (—00 < ^ < s) and notice the property of the Pred- 
holm determinant det(l — Q1Q2) = det(l — Q2Qi)- We easily find that (52Qi(Ci)C2) — 
Xo(6)^x(6,6)Xo(6), where i?x(6,6) is given in (2.10). 

Because an approximation is involved in the derivation, the expression (2.25) is most 
likely not exact when n > 2. A validity of employing this approximation is that, for the 
narrow wedge case, the generating function becomes the multi-point distribution function 
for the Airya process in the long-time limit which is expected from the universality. In 
this paper we will see that the situation is similar for the case of the half Brownian motion 
initial condition. In the long-time limit, the multi-point generating function is equivalent 
to the multi-point height distribution function, 

lim G^,{{s}n, {X}n) = lim Prob (ht{Xj) < sj, j = 1, • • • , n) . (2.33) 
From (2.25), we find 

jim Gl {{s}n, {X}„) = det (1 - P,/Ci2P.) 



M 




M\ 11 \ ^ / ^^^^^"^ ^^'^ (^12(^n,, e.-; (2.34) 

M=0 ■ 1=1 \ni=l "^""^ / 

where Xs{C) is defined below (2.32) and 

lCi2{Xi, {1; X2, .^2) = /C2(Xi, ^1; X2, ^2) 

+ Ai(e2) (^e-^+^^«^ - dye-''-y Ai{i2 + y)^ , (2.35) 

[-J-00 ' "'Mij + w)Ai{^k + J <k. 

The Fredholm determinant with this kernel has appeared in the PNG model and TASEP [42, 
43, 44]. Note that for the one-point case n = 1, it reduces to (2.22). 
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As pointed out in the Conjecture 9 in [38], it is plausible that in the long-time limit, the 
Fredholm determinant det (1 — PsK,i2Ps) describes the multi-point distribution considering 
the exact result of the one-point case. Thus (2.34) strongly suggests that the factorization 
approximation becomes exact in the long time limit in both narrow wedge and half Brow- 
nian motion initial conditions. In Sec. 6, we briefly discuss the derivation of these two 
equations (2.25) and (2.34). 



3 KPZ equation and (^-Bose gas 
3.1 Feynman-Kac formula 

In subsection 2.1, we mentioned that the KPZ equation (1.1) is not well-defined as it is, 
and defined its sohition by (2.1) and (2.2). Here we explain this using a limit of a modified 
version of the KPZ equation. 



dhf,{x,t) A / dh^{x,t) 



f dKix,t) Y , d'K{x,t) , ^ U^-^Vntn^ fo,^ 



dt 2 



This looks almost the same as the original (1.1), but there are two differences. First rj^ix, t) 
is still a Gaussian noise but its covariance in x-direction is now smeared as 

{V.ix,t)7]^{x',t')) = Cix - x')Sit - t'), (3.2) 

where C^ix) is written as C^{x) = k,C{kx) and C{x) is a smooth, even and positive function 
such that J^^C{x) = 1. Second, the constant velocity term — (^\^/D/2uj Ck(0)/2 is 

added. Note that the above properties of imply lim^^ooCK{x) = S{x) so that the 
noise 77^ tends to the white noise in the original KPZ equation. An apparent problem is 
that the additional term in (3.1) diverges as 5(0) but in fact this term is necessary for 
considering a meaningful limit of the equation (3.1). To see this let us apply the Cole-Hopf 
transformation, 

A 



Z^(x, t) = exp y^K(x, t)j , (3.3) 
to (3.1). By Ito's formula (3.1) becomes 

dt ~ dx^ + ^^-^^ 

which is a well-defined stochastic differential equation of Ito-type. Clearly the solution 
Z^ix^t) to this equation converges to that Z,y^x,D{x,t) of (2.2) as k — >■ 00 [20]. Hence 
the Cole-Hopf solution (2.1) with (2.2) can be interpreted as a hmit of the solution to the 
equation (3.1): 

^i/,A,^(a;,^) = lim ^K(a;,t). (3.5) 
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Note that if the constant velocity term is absent in (3.1), there appears an additional term 
in (3.4), or the equation (3.4) should be understood as representing a stochastic differential 
equation of Stratonovich type, which is not well-defined in the limit /t — )> oo. 
The solution of (3.4) can be represented as the Feynman-Kac formula [46], 



ZK{x,t) = Ero exp 



2iy 



T]^ {h{2us),t — s) ds 



Z,(6(t),0) e" 



1 ( X-/D 



Cn{0)t 



, (3.6) 



where represents the averaging over the standard Brownian motion 6(s), < s < t 
with 6(0) = X and the initial condition is 



Z,{x,t = 0) 



0, 



X <0, 



,aB{x)^ X > 0. 



(3.7) 



Here B{x) represents the one-dimensional Brownian motion with 5(0) ~ 0. Some readers 
may find the following path-integral expression more intuitive. 



Z,{x,t)= r dy r^'^ ^ D[x{T)]exp\-S[x{T)]+aB{y)-l{^\ C,{0)t \ , 

Jo Jx{0)=y 



2v 



where the action 5'[x(t)] is given by 



SWx)l = /'..(i-( 



1 / dx{T) 



dT 



2v 



(3.8) 



(3.9) 



We find that Z^x,D{x,t) defined in (2.1) is written in terms of Zi ^ i{x,t) with specific 
parameters u = 1/2, X = D = las follows. First using the Feynman-Kac formula (3.6), 
one has 



= lim Ej. I exp 



Av^ f^^ f..^ . t \ 1 ( Av^ \ ^ ^ 
r]^[b{2iys), s] ds - - — — C^{0) — 



2v 



2v 



2v 



2v 



Z(°)(6(t),a)j , 
(3.10) 



where we rewrite ^^(a;, 0) as Z'^'^\x,a) in order to represent explicitly the dependence on 
a and the independence of k. From the properties of written below (3.2) we find the 
scaling relations of and /^^(a;, t), 



C^{x) = aC^/a{ax), 
r]^{x,t) = {aby/'^r]^,/a{ax,bt), 



(3.11) 
(3.12) 
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where a, 6 e R and the equahty (3.12) holds in the sense of distribution. By use of (3.12) 
with a — 1,6= l/2i^, (3.10) becomes 



hm exp a / r]^ {b{s),t - s)ds - —C^{0)t 
\ I Jo ^ . 

Next one remembers the scahng property of the Brownian motion, 



(3.13) 



(3.14) 



where the equahty holds again in the sense of distribution and wc put the initial position 
X of b{s) explicitly. By (3.12) and (3.14) with a — a'^,b — a'^, Eq. (3.13) can be rewritten 
as 



lim ^a'^x cxp 



J Vn/a^ {b{s), a t- s)ds- YC«/a2(0)t 



Z^'>\a-''b{aH),a) 



At last noticing from (3.7) and (3.14) that 

Z^'^^a-^'x^a) = Z(°)(x,a = 1), 

we find (3.15) is nothing but Zi ^ ^(a^x, aH). 
Thus we have established 



(3.15) 



(3.16) 



(3.17) 



In what follows we restrict our discussions to hi i i{x, t) and Zi^^^^{x, t) (hereafter we omit 
the indices |, 1, 1). We remark that, for this special parameter values, (2.3) reads 



7t 



, x = 2^^X. 



(3.18) 



3.2 The (5-function Bose gas 

Next we consider the replica partition function {Z^{x^t)) {N = 0, 1,2, ■■■). First we 
perform the Gaussian average over riK{x,t) using the path integral representation (3.8) 
and (3.2), to find 



(3.19) 



The right hand side of this equation includes the self-interaction term '^f=k=i ^ni^j —Xk) — 
NCk{0), which becomes divergent in the limit k — > oo. This term, however, cancels out 
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the last term in (3.8), and thus we can take the hmit k — >■ oo. We obtain 

(Z^(x,t))= hm(Zf(a;,t)) 



^ poo P 
i=l -^0 Jx, 



Xj{t)=X 



(0)=% 



D[xj{T)] exp 



N 



dT 



3+k=\ 



N 



X ( exp J^5(|/fc) 



(3.20) 



where (■ • • ) in the last factor indicates the remaining average over the Brownian motion 

B{y). 

The right hand side of this equation represents the imaginary-time dynamics of the 
5-function Bose gas with attractive interaction with the Hamiltonian Hn, 



N 



H 



N 



2^9x2 



in terms of which the replica partition function {Z^{x, t)) can be written as 



{Z^{x,t)) = {x\e-''^'\^) 



—Hivti 



(3.21) 



(3.22) 



Here {x\ represents the state with all N particles being at the position x and the |$) the 
initial state of the 5-function Bose gas. One can perform the average over the Brownian 
motion B{y) and the dependence of |$) on xi, . . . , xn can be explicitly calculated as 



■ PeSN \ \k=i / / 



Pi 

N 

En 

PeSn j=i 

N 



( 



\ 



dzi 



V'27r(xp(j) - xp(j_i)) 



e {xp^j) - a:p(j_i)) 



(3.23) 



Here Sn denotes the set of permutations of order A^, Q{x) is the step function and we set 
xp{o) = 0. Since we are considering a Boson system, the above function is taken to be 
symmetric in Xi, ■ ■ ■ ,Xn- 

The eigenvalues and eigenfunctions of the (5-Bosc gas can be constructed by using the 
Bethe ansatz [27, 28, 29, 31, 32]. Let 1^^) and be the eigenstate and its eigenvalue of 



H 



N, 



(3.24) 
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By the Bethe ansatz, they are given as 



(xi, • • • , xjvl^z) = C'^ sgnP Yl i^PU) - ^P{k) + 'isgn{xj - Xk)) exp ( / ^ zp(i)Xi 

PeSN l<j<k<N \ 1=1 

(3.25) 

where Cz is the normahzation constant, for which a formula is given in (3.27) below. 

For the (5-Bose gas with attractive interaction, the quasimomenta Zj {1 < j < N) which 
label the state, are in general complex numbers, zj (1 < j < A^) are divided into M groups 
where 1 < M < N. The ath group consists of n„ quasimomenta ZjS which share the 

common real part Qa- Note that X^aLi = ^- The quasimomenta in each group line 
up with regular intervals with unit length along the imaginary direction. Using qa and 
n-a < 0£ < M), we represent zj {1 < j < N) as 



a-l 



^3 



?a - ^ ("-a + 1 - 2rc,) , for j ^^n^ + r^, (3.26) 



where 1 < a < M and 1 < < no,. The normahzation constant Cz, which is taken to be 
a positive real number, and the eigenvalue are given by [28] 



l<j<k<N ' ' 



N ^ M M 



j=l a=l a=l 

There is a problem of completeness of these states. This has not been resolved rigorously 
yet but the fact that one can recover the height distribution of the KPZ equation for the 
narrow wedge initial condition is a strong affirmative evidence. Here we proceed assuming 
its vahdity and will see the consistency with the computations based on the ASEP. This 
provides a further evidence that the above Bethe states are in fact complete. 



4 Generating function 

In this section, we give a derivation of the formula (2.6) in Theorem 1. The replica partition 
function {Z^{x,t)) (3.22) can be written as 

/oo /»oo 
dyi--- dyN{x\e-"^'\yi,--- ,yN){yi,--- ,yN\^). (4.1) 
■oo J —oo 
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Expanding the propagator {x\e • • • , ^at) by the Bethe eigenstates of the (5-Bose gas 

(3.25), we have 

M=l j=\'>-°° \-' a=\ na=\) 

X e-^^*(a;|^,)(^,|yi,--- ,yjv)(yi,-- - ,l/iv|$). (4.2) 
Here we want to perform the integrations over (1 < j < A'"), 

^ poo 

n/ ,2/iv)(|/i,--- ,l/iv|$) (4.3) 

using (3.23) and (3.25). But this is not allowed for the moment because the integra- 
tions over Qa, (1 < a < M) must be performed before those over yj, {1 < j < N). 
To see this explicitly, one uses (3.23) and notices the symmetry of the eigenf unction, 
{'^z\yp(i), ■ ■ ■ ,yp{N)) — {^z\yi, ■ ■ ■ ^Vn), to find that the right hand side of (4.2) is repre- 
sented as 

An /-oo / poo J ~ \ 

M=l j=l"^%-i a=l V-^-^ " na=l/ 

TV 

X e-^^*(a;|^,)(^,||/i,-- - ,|/^)ne^^'^"'''^'^'^- (4-4) 

It is clear that the integrand on the right hand side of this equation is not integrable 
on Vj {1 < j < M) due to the factor exp((2A'" — 2j -\- l)yj) while it is integrable on 
?a (1 < a < M) thanks to the factor e-""^'*/^ in e'^'K 

4.1 Deformation of contours 

To exchange the order of integrations, we notice 

Lemma 4 In (4.4), we can deform the contour of qa {1 < ck < M) to R — ic where c is an 
arbitrary real constant. 

Proof In (4.4), we change the variables from to (1 < J < M) such that 

Ma = ga - ga+l, (1 < < M - 1), UM = qM- (4.5) 

The singularity in Ua of the integrand in (4.4) comes only from the factor \Cz\^ ■ Note that 
it depends only on mi, • • • , um-i and is independent of um- Thus there are no poles of um 
in (4.4), and we can deform the contour of um to R — ic. Noting that 

Yi f ^ = fr / ^, (4.6) 

H loo 27r 27r £1 J^_,, 27r ' ^ ' 

we find the statement holds. 
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Using Lemma 4, we deform the contours of g^, (1 < a < M) to R — ic where c is a 
positive constant so large that the factor exp (^—i X^^^ ^p(j)yj + I Z^jLi(2-^ — 2j + l)yj^ 
in the integrand of (4.4) converge when yj > goes to infinity. Here Zj = qa(j) + 
t (^a(j) + 1 — 2r(j)) but one has to be careful about the notation. This is a usual complex 
conjugate of Zj before lemma 4, i.e., when q^Q^ is real. But now it is not because after the 
deformation of the contour, q^^j) has an imaginary part. 

After this deformation, we can perform the integrations of yj {1 < j < N) before those 
oi qa {1 < a < M). Eq. (4.2) can now be expressed as 

N M / „ , oo\ 

M=l ■ a=l \'^K-ic n„=lj 

Here {x\^z) is given by (3.25) with xi — ■ ■ ■ — xn — x and is computed as 

N r-OO 

(^,1$) = JJ / dyj{m,\yi,--- ,|/iv)(2/i,--- ,2/iv|$) 
= J^sgnPn/ dz/.e-<^^«+t(--^+^))- n 



Pe5jv i=l ^'-1 l<j<k<N 

^ -1 

^ivic.J^sgnP n (4o-)-4w + ^)n _,(,. + yT7v2- 

(4.8) 

In (4.7), we take the imaginary part c of q^ in such a way that we can perform the 

integrations of yj (1 < j < N) in (4.8), i.e. Re(--i(zp(^-) H h -2p(Ar_,+i)) + /^/2) < for 

any / (1 < / < A^). For example, if we fix c such that c > N/2 + maxQ,nQ,/2, the above 
condition is satisfied. 

4.2 Combinatorial identities 

For further analysis of the integrand in (4.7), we need two combinatorial identities for 
{x\^z) and The first one is for (xl^^)- One has 

j2 sgnp n (^^0) - + ^)) = n - ^'^) (^-9) 

Pe5jv l<j<k<N l<j<k<N 

for any complex variables Wj {1 < j < N) and f{j,k). This identity was derived as 
Lemma 1 in [33]. 

The next one for the term is 
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Lemma 5 For any complex numbers Wj {1 < j < N) and a, 

^ 1 
V sgnP TT (wp( •) - wp(fc) + a) TT ■ ■ ■ ^— 

^ 1 

l<j<k<N m=l ^ ' 

Proof A similar identity appears in the context of ASEP with the step BernouUi initial 
condition and has been proved in section 111 of [19]. Here we follow the same strategy using 
mathematical induction. Let us call the left hand side and the right hand side of (4.10) as 
ipj^ and respectively, i.e., 

'0jv(tyi, • • • , wn) := ^ SgnP {'^P{j) - '^P{k) + a) 

PeSN l<j<k<N 

^ 1 

X TT (4-11) 

wp{N) H h wp(N-m+i) + m^a/2 

^ 1 

0jv(wi, • • • , wn) := TT ^^J' ~ '^''^ TT T~lo- ^^--^^^ 

l<i<fe<iV ^=1"^"^ + "/^ 

We want to show tjj^ = (f)^. We easily see that it holds for the case = 1. Let us assume 
that it holds for A^ — 1 case. In (4.11), we first sum over all permutations with P(l) = / 
fixed and then sum over I. Noticing Yli<k<Ni'^i ~ 'Wp(k) + a) = nfc(^i) {''^i ~Wk + a), we see 

iPn{wi, ■ ■ ■ ,Wn) 

1 ^ 

^(-1)^+^ TT ~ ^-J' + o,)'4'N-i{wi, • • • , Wi-i, Wi+i, ■■■ , Wn) 



N 



^(-1)'+^ TT ^"^^ ~ '^j + a)(j)N-i{wi, ■■■ , wi^i, wi+i, ■■■ , Wn), 



(4.13) 



where we used the assumption of the mathematical induction, ■0iv-i — <f>N-i, in the second 
equality. The identity ipN — 4>n is now equivalent to 

Ei 1^^+lTT^ , ^ </>Af-i(w^i, ■ ■ ■ , wz+i, • • • , wjv) A a 
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Prom the definition (4.12) of 0Ar, we get 

4>N-l{wi, ■■■ , Wi-i, Wi+i, ■■■ , Wn) 



Wi + a/2 



U^Jiiwj - wi) ■ Ylf=i+ii'Wi - '^j) 



0Ar(wi, • • • ,Wn) 



wi + a/2 



wi, • • • ,wn)- (4.15) 



Substituting this equation into (4.14), one sees that it is now enough to show 

E(--<'/2)n 1 + ^ =E'". + ^- 

1=1 ji^l) ^ j=l 

This is proved as follows. One notices that the left hand side can be represented by a 
contour integral, 



(4.16) 



/ dz(z + a/2)f[(l + ^-), 
2maJc„ |j^V ^-^,7 



(4.17) 



where Cr is a contour enclosing the origin anticlockwise with radius R taken to be so large 
that the contour surrounds all the poles wj {1 < j < N) in the integrand. The contour 
integration can be performed as follows. Expanding the product in the integrand, one finds 



l<l<m<N 
a 



{Z - Wl){z - Wni) 



a\ 
2) 



■^■^ z — Wi V 



(4.18) 



■ 1 - -'J 

Further expanding the second and third terms in the integrand in 1/ z, we see 



z-wi 



z + 



2/ 



a + 



a{wi + a/2) 



+ 0(^1, 



(a 

r- Z H 

{Z - Wl)[z - Wm) V 2 



+ 



(4.19) 



The higher terms in (4.18) are of order 0(1/2;^) and thus do not contribute to the contour 
integral. Hence (4.17) is calculated as 



/ dz{z + a/2) f[(l + 
27naJc^ z-WjJ 



2ma 

N 



N 



Cr 



1=1 ^ l<l<m<N ^ 



N 



1=1 



l<m 



2 ' 



(4.20) 



thus we obtain (4.16). 



19 



Using the identities (4.9) and (4.10) to (3.25) and (4.8) respectively, we get 

(x|*,)=7V!C, n {zj-Zk)e'^"=^''-, (4.21) 

l<j<k<N 

^ 1 

= r^Tvic, n (^l-^^)n 7^-77^- (4.22) 

l<j<k<N 1=1 I ^ ' 

Thus using these relations and (3.27), we find that the factor {x\^z){^z\^) in (4.7) becomes 



[x 



a=l l<Kfe<ivl^^ ^1 Z=l^^' 



Here we used the fact that z* — zl is the complex conjugate of Zj — z^ although z* is not 
that of Zj as mentioned below Lemma 4. The common imaginary part c of go,(j) and qa(k) 
cancels out by the subtraction. 

We want to rewrite this equation in terms of and Ua in (3.26). For the last factor of 
this equation, we easily find 

From (3.27) and (4.21), we know the remaining factors in (4.23) is represented by |(0|\I'2)p. 
For this quantity, the following result was obtained in Appendix B in [28], 

• ) "I r I I 



i(oiv^.)r=iV!n^ n 



a=l l<j<k<N ' ■' ' 



N\ T-r |ga-g^-|(^a-ri^)r ^.x 



From (4.23)-(4.25), we get 



We can further deform (4.26) to an expression in terms of a determinant by using the 
Cauchy's determinant formula, 

nf</3(«a -^/3) . M(M-i) / 1 \ 

=pvl — r\ = (-1) ' det , (4.27) 
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and a few basic properties of determinant. We find 



\ M M ria 

— ) rr e'^""^"" TT — ^- 

Hj + nk + 2i {Qj -qk)J - i^^^ ^J^ iQa + i (n„ - 2r ) 

1 

[ria - 2r) 

= 2^N\ TT / duja] det e'"^*^-"^^'"^^'^''^-^^*^'^^ TT w r 



2^^A^! IT ( / ^'^^ ) det (e-'^^K+"'=+2^(9.-9fc))) ^^^^ JT e^""^"^ JT 

a=l / ^' a=l r=l + 2(^« 



(4.28) 

where in the last equahty, we used a simple fact 

det(a5^+'"=) = det{{ajakp). (4.29) 

4.3 Fredholm determinant representation of the generating func- 
tion 

Prom (3.28), (4.7), and (4.28), we obtain an expression of {Z^{x,t)) in terms of the deter- 
minant, 

{Z^{x,t))e^^— 

M=l a=l \na=l"^"* / 



X det 



2 \ M 

girij qjX+rij ^ -Uj {ojj +uik)-2igj {u}j -uj^. ) 



n"ii% + K%-2r) 



M=l a=l V"^" na=l/ 



X det ' ' ^ 



n;iii?+i(n,-2r) 



(4.30) 



The imaginary part c of the contour of q is discussed below (4.7) and wc find it is sufficient 
to satisfy the condition c > N/2 + maxana/2. Note that in (4.30), the contour satisfying 
this condition passes below the poles q — j — na/2 (j = 0, 1, • • • , [na/2]) where [a] is the 
largest integer which is smaller than a and that we can deform the contour as long as it 
satisfies this property. Thus in (4.30), we can relax the condition to c > maxQ,na/2. 



21 



Substituting this equation into (2.5), we eventually get the Predholm determinant rep- 
resentation of the generating function (2.6). Remembering (3.18), we have 



Ar=lM=l ■ a=l V"^" na=l/ 



X det 



^3 \ M 



^;ll(^?+|(r^,•-2r)) 



j,k=i 



M=0 a=l " 



X det y (-1)"-' / —- FF^^- n TV 

[ti -^K-^c^^ n:=i(^9+K^-2r)) 

with Cn > n/2. Shifting the variable q to q + iX/^t, we get 

^ XI A/fi n 



j,k=l 

(4.31) 



M! 

M=0 ik=l 



^ JM-ic„ TT nr=i(«?- :^ + i(^-2r)) 



with Cn > X/jt + n/2. Applying the relation 

1 ^(iq---^) 

TT t = V ^* y u 33) 

l},iq-Xh, + \{n-2r) r (i^ - | + f ) 

to (4.32), we finally obtain (2.6). 



5 Another expression of the kernel 

In this section we give a derivation of the expression (2.10) in Proposition 2. We use the 
following two relations. 

Lemma 6 
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(a) We set a e IR and m, n > 0. When Im q < —n/2 + a, we have 

r (^g + g - f ) 
r (ig + a + f ) 

where 



e 3 



y^dyAi? (^l/,^,ig + a)e--^ (5.1) 



A^K,e)^^|^..e--/3^£^ (5.2) 



Fp defined in (2.12) below. 
(b) For u,v,x gM. and w > 0, we have 



23 



1 /"°° 

^TJ" ,/-oo 

^AF (2-l{v + x),2^w,uj Air (2-^v - x),2^w,uj , (5.3) 



where Ai (a, 6, c) and Air(a, 6, c) are defined by (2.11) and (2.12) respectively. 
Proof 

(a) The right hand side of (5.1) is written as 



rH2^ + ^? + «) 



Here we used the fact that the contour of z can be deformed from T2m(ia-q) to M + imn 
since the imaginary part of the poles 2m{ia — q + ir), r = 0,1,2, ■■ ■ of the integrand is 
larger than mn when the condition Img<a — n/2is satisfied. In this equation, we change 
the variable z on the right hand side to y2 = — iinn and get 



271 ' T{iq + a + l-i^) 



oo 



S(y r(^g + a-f + ^f^) _ r(^g + a-f) 

r (zg + a + 2 - zf^) r[zq + a + ^) 



(b) The left hand side of (5.3) reads 



^ ^^dp[ dze^^^^^^^^^^^ T{nr{^+p) + a) ^^^^^ 



Jt_^^,^ V {-iw{z - p) u)' 



By applying the change of variables p — {zi — Z2)/2'^/^ and z — {zi + Z2)/2'^/^, we obtain 
the desired expression. 
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Using Lemma 6, we can obtain (2.10). Applying (a) of the lemma with a — —X/^t and 
m = 7t/22/3 to (2.7), one has 

Kx{u,,u,) = f;(-l)"-^ / ^ r dy.kil ^,iq - ^) e-2^^(-^--^) 



At this point, the contour of q can be replaced to R since since the contour of z in the 
definition of ^r(a, fc, c) (5.2) pass below the singularity of the Gamma function. Changing 
the variables g, and top = 2^/^7ig, = 2j^^ujj and y = yi/2^/^-p^/2^/^-(^j-+^fe)/2, 
we see Kx{ujj,uJk)doJj = Kx{Cj,Ck)d^j, where 

^TI" J-OO ^-OO \ 23 237i 237t J 

oo 

X ^(-l)"-ie-^*"(^-2'). (5.8) 

n=l 

Thus we can apply (2) of Lemma 6 to this equation with v = 2^/^y + {^j + ^jfc)/2^/^, 
w = 1/(2^37^), u = -X/-ft, X = {^j - Cifc)/2^/^ and arrive at (2.10). 

6 Multi-point distribution 

6.1 Generating function 

We consider the n-point generating function defined by (2.24) and derive the result (2.25). 
Eq. (2.24) can be expanded in terms of the replica partition function {Z{xi-^,t) ■ ■ ■ Z{xi^,t)), 

(— Ij e 24 



G,,{{s}n,{X}n) = l + J2 m 

N=l 

X 



e-^^^.-^(^^-^f\Z{2j^Xi„t) ■ ■■Z{2j^Xi,,t)). (6.1) 

As we discussed in Sec. 3 and 4, we can express the replica partition function in terms of 
the Bethe states (3.25)-(3.28) of the 5-Bose gas. Specifically we can the generalize (4.7) to 
the case where x's are distinct and we get 

{Z(xi„t)---Z(xi^,t)) 

N M / „ , oc\ 

M=l ■ a=l \'^K-ic no,=l/ 
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Thus the generating function is written as 



G.. {x}n) = 1 + E n / . It ^ ^^f^^ "--^^^^ 

AT 1 * „ 1 \ t/ M. — iC „ 1 / 



Ar=l a=l V^K-JC -■■ „„=iy 

where 



A, = e-^^*(0|^',)(*,|$), (6.4) 

(01*) ■ ^ ■ ^ 

Note that in two factors Az and Bz, only Az depends on the initial condition and we can 
readily use (3.28) and (4.28) for Az- On the other hand, the factor Bz, which is irrelevant 
to specific initial conditions, has already appeared in the narrow wedge initial condition 
and the following approximation has been proposed in [33, 34]. 

For Bz, we have to estimate the equation with the following type, 

n N 

{xi„--- ,xi,\^z)lle'^'^. (6.6) 

h,---,lN=l 3=1 

Here gj {j = I,-- - ,n) are variables which depend on j. Although the eigenfunction 
{xi-^, ■ ■ ■ ,xi^\'i/z) is given in (3.25), the following compact form given in [28] is convenient 
for our current discussion, 

(M ^ M \ 

i^Qa y^'^Yl \yc-yc'\\, 
a=i cena(P) a=i c,c'ena(P) J 

(6.7) 

where Ap, fla{P) and the prime in the summation are defined as follows: The amplitude 
Ap is represented as 

Ap{yi, ■■■ ,yN)^ sgnP • Cz Yl{k\) (zp(j) - zp(^k) + isgn{yj - yu)). (6.8) 

a=l fe=l l<j<k<N 

The cluster Qa(-P) is defined by ^a{P) = {a\a{P{a)) = /3} by use of the cluster counting 



function a : [1, ■ ■ • , A^] [1, ■ ■ ■ , M] such that a{a) = /3, for ^'^^^ nj < a < ^J^^ 
The prime means the summation over permutations which keep the order inside each 
cluster, i.e., the permutations satisfying the condition j < k ^ P{j) < for j and k 
such that a{P{j)) = a{P{k)) (1 <j,k<N). 
Thus (6.6) is expressed as 

n N n 

Y {xi„--- ,x,,\^z)ll'''^ = E Y^Apixi,,--- ,xi,)ef(i^^-'^\ (6.9) 

'i,---i^Ar=l J=l h,--- ,In=1 P&Sn 
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where (f){{l}N, P) represents the phase 

M ^ M N 

(t>{{i}N,P)^i'^qa ^ Vc-^Y^ \yc-yc'\ + Ygi, (6.10) 

The difficulty in (6.9) is that the summations over Zi, • • • ,In and P are coupled. In [33, 34], 
the authors proposed an approximation that the two summations factorize, 

n n 

Y: j:'Ap(y,„...,y,Je^(«.,P)^ ^ e^«'>-) ^^'^^(y,,, • • • , (6.11) 

where 0({Z}iv) = 0({^}iv, -P = (1, 2, • • • , N)). Under this assumption, we can perform the 
two summations using the relations 

Y'MVhr-- ,yir.) = m.), (6.12) 

P£Sn 

which comes from the identity (4.9), and 

n M / n \ 

^<t'i{i}N) ^Y^^-lY:lk=i\=^j-^k\99jd9k ('^(^m+iqc.xA ^ (6.13) 

h,---,lN=l a=l \l=l / 

which is obtained as (3.13) in [34]. From (6.11)-(6.13) and setting yj — 2^^Xj, gj — sj—Xj, 
we obtain the approximated form of (6.5), 

a=l \/=l / 

Substituting the equation above into (6.3) and following the same procedure in Sec.4.3, we 
have 



{xu = 1 + E ]^ n / E W-Ht 



1 



(ria + np) + i{qa- qp) 



M n, 

X 



nn V^g, + i(n„-2r), 

M / 2 " ^ 

a=l \ Z=l 

(6.15) 

Here the index jj indicates that we use the approximation (6.14). We find that the equation 
above has a similar structure to the corresponding one in the narrow wedge case in [34]. 
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All the differences are just two points: The contour of q includes the imaginary term —ic 
and the factor 11"=! ( ig -2r) ) added. Thus we can just follow Sec. 4 in [34] 

if we replace Eqs. (4.14) and (4.19) in [34] by (1) and (2) in Lemma 6 respectively. As a 
result, we find the Predholm determinant expression of G^^({s}„, {X}n): 

GU{s}n, {X}„) = det(l - L), (6.16) 



7t 

where 



(-|^ n n n \ 

~2 ^ ~ ^'^l^^^-^^fc - 7t XI Xpsj - {dz - dz') Xjdsj j 
j,k=l j=l j=l / 

X J duAi^ (u + z,^,0^ Air (^u + ^{{u- s}„). 

(6.17) 

Here Xs{x) is given below (2.32), Ai'"(a, b, c) and Air(a, 6, c) are defined by (2.11) and (2.12) 
respectively and for $({a;}„), see (2.27). This expression of the kernel corresponds to 
Eq. (4.25) in [34] for the narrow wedge initial condition. The only difference is that 
the product of the ordinary Airy functions Ai(u + z)Ai(u + z') is replaced by Ai^(u + 
z, l/7i, 0)Air(ix + z', \hu 0) in (6.17). 

We can further deform the expression (6.17) following the discussion in Sec. 4.2 in [34]. 
The equation corresponding to Eq. (4.38) in [34] becomes 

/oo 
rfMi • • • dii„(Ai^,,i/^^,_^„/^jLoe(^°-^i)^|iii)(iii|e(^i-^^)^|ii2) x • • • 

X (M„_i|e^^"-i-^-)^|Mn)K|e(^"-^°)^Lo|Air,.,i/^,,-Xo/7t)'^({w - s}n), 

(6.18) 

with — 0. Here B. is the Airy Hamiltonian given below (2.26), $({x}„) and Lq are 
defined by (2.27) and (2.28) respectively. The states |Ai^_j,_^) and |Air,z,6,c) are defined by 
the relations (M|Ai^^^) = Ai'"(2; + m, 6, c) and (w|Air,z,6,c) — Air(-z + 6, c) respectively. In 
this kernel, we put L = Q2Q1, where 

Q2{z',m) = (Ai^,_i/^^,_;,„/^jLoe-^i^|i.i), (6.19) 

/oo 
du2 ■ ■■dun{ui\e^^'-^^^"\u2) X • • • X (M„_i|e(^"-i-^")^|M„) 
-00 

X (K„|e^"^Lo|Air,^,i/^,,-Xo/7*)*({^* - (6-20) 

Noticing the relation of the Predholm determinant, det(l — Q2Q1) = dct(l —Q1Q2) and the 
biorthogonahty (A. 2) of the deformed Airy function, we get G^^({s}„, {X}n) — det(l — Q) 
where 

Q{x,y) := QiQ2{x,y) 

/oo 
du2 ■ • •dxi„(iii|e(^^-^^)^|ii2) • • • K|e^"^Lge-^^^|y)$({ii - s}„). (6.21) 
■00 
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The factor {x\e^-^^~-^'^'>^\y) is represented in terms of the deformed Airy functions as (2.29), 
which will be shown in Appendix. Using (A. 2) and (ii) in the Appendix, we find Lq = 
Lq and Lqc^^^^ = e^^^^Li, and thus we finally obtain (2.26). The second last factor 
{x\e''^ Lj\y) in (2.26) is written as 

(x|e*^L,.|y) = r dze-'^kir (x + -, Z^k±l\ Af (y + ^, 1, . (6.22) 

Jo \ It It J \ It It J 

This readily follows for i < from (A. 4) but one has to be careful when t > 0. The 
integrand is divergent as ^ ^ — oo and hence {x\e^^\y) in (2.29) is not well-defined. But 
in {x\e^^ Lj\y), Lj projects the range of integration to the positive side only and (6.22) is 
still valid. 



6.2 Long-time limit 

Here we consider the long-time limit of the generating function G^^({s}„, {X}„) (2.26). 
We first introduce the limiting version of Lj{x,y) (2.28), 

POO 

Cj{x,y) :— liui Lj{x,y) = / dwAr {x + w, Xj)Ai^ {y + w, Xj), (6.23) 
where Ai^{x + w, Xj) is the limit of the Gamma-deformed Airy functions (2.11) and (2.12), 

1 / I —v\ 1 f ^izx+if- 

Ar(x,y) lim — Ai^ x, — , — = — / dz (6.24) 

7t-^oo 7t y It It J 27r Jy_^^ iz-y 

Ai+{x,y) := lim 7tAir ( x, -, — ^ = — / dze''''+'^ {-iz - y). (6.25) 
Tt^oo V It It J 27r 

Here the contour V_iy in (6.24) is given below (2.12). Noticing 

n 

lim ^{{u - s}n) = 1 - TT (1 - Q{sj - Uj)) , (6.26) 

where Q{x) = 1 (a; > 0), (x < 0), we find the long-time limit of (2.26) becomes 

lim Gl{{sU {X}n) = det (l - £i + P.^e^^-^^)^ • • • P.„_,e(^"— ^")^P,„e(^-^^)^£i) . 

(6.27) 

Here P^ represents projection operator onto (— oo, s). To this equation, we apply the result 
of Appendix in [34] with = 5ijK replaced by K^^ — 5ijCj. We eventually obtain the 
Predholm determinant with matrix kernel M. 

lim GU{s]n, {X]n) = det (1 - M) , (6.28) 
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where 

Here Xs(a;) is given in (2.34) below. 

We find that this matrix kernel is equivalent to ICi2{Xj, ^i, X^, ^2)- By taking the limit 
7t — )■ 00 in (6.22), we easily find the factor {^i\e^^i~^''^^ £k\^k) in the case j >k oi (6.29) 
is expressed as 

/•oo 

= / (i^e-(^^-^'=)^Ai-(ei + z, X,-)Ai+(6 + z, Xk) 
Jo 

J_oo J-00 («^1 - ^j) («^1 + - + ^fe) 

Here we used (6.24) and (6.25) in the second equality. Noticing 

+ Xk ^ 1 1 g^V 

{iwi - Xj){iwi + iw2 - Xj + Xk) itih - Xj iwi + iw2 - Xj + Xk' 
we finally find the factor can be represented as the Airy function. For j > k, we obtain 

POO 

(6|e(^^-^'=)^A|6) = / d^/;e-"'(^^-^'=)Ai(6 + w)Ai{^2 + w) 

Jo 

+ Ai(6) ^e-^'+^^^i - d«;e-^^"'Ai(ei + w)^ . (6.32) 
Similarly for the case j < A; in (6.29), we get 

(6|e(^^-^'=)^A|6) = - f dwe-"'(^^-^'=)Ai(6 + w)Ai{i2 + w) 

J —00 

+ Ai(6) (^e-^l+^^-^i - d«;e-^^^Ai(ei + w)^ . (6.33) 
Hence we finally find Mjk{^i, ^2) = Kii2{Xj, ^1; X^, ^2)- 

7 Directed polymer interpretation 

Our result can also be rephrased for the free energy distribution of the 1+1 dimensional 
directed polymer in random media. First let us rewrite the path integral expression of 
Zy,\,D{x,t) (3.8) as 



Zis,jA^, t)^ dy D[x{r)] exp [-P(H[x] + iiB(y))] , (7.1) 

^0 Jx{0)=y 
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where 




VDr]{x{r),T) 



(7.2) 



The parameters /3, 7 and are defined as 




(7.3) 



in terms of those of the KPZ equation (1.1). 

This expression (7.1) can be interpreted as the partition function of a directed polymer 
in random media with inverse temperature /3. On the two dimensional {x,t) (x G R and 
r > 0) plane the configuration of a polymer is represented as a function x{t) {0 < t <t). 
Here it is assumed that x{t) is a single- valued function, which means that the polymer 
is "directed" for r coordinate. For a configuration x{t) (O < r < t) the bulk energy 
H[x] (7.2) is assigned. The first term of H[x] represents the elastic energy of the polymer 
whose strength is determined by 7, and the second term represents a random potential 
energy with D adjusting its strength. 

The half Brownian motion initial condition (1.4) in the KPZ equation corresponds to 
the following boundary condition in the language of the random directed polymer. At 
the boundary r = 0, the position of the polymer x{0) can take any positive value y > 
and, depending on the position y, the random boundary energy fJ.B{y) is assigned to the 
polymer. Note that the other edge x{t) is pinned at x. As a statistical mechanical model, 
H in (7.1) is an independent parameter, but when translated from the KPZ equation one 
has to keep in mind that it is written as ^ = ^/P^-yD. 

We want to discuss the low temperature behavior of the directed polymer for a fixed 
t. When the temperature is sufficiently low, the configuration which gives the minimum 
energy is dominant for the free energy. When the position x of the edge x{t) of the polymer 
is positive, the other position of the edge y tends to be close to x because otherwise the 
polymer obtains excess bulk elastic energy. In this case, it is dominated by the boundary 
energy B{y). Thus in the positive region a; > 0, the free energy fluctuation is described by 
the Gaussian. On the other hand, when x < 0, the edge x(0) tends to be pinned at the 
origin. In this case, the free energy fluctuation is determined only by the bulk energy (note 
that 5(1/ = 0) = 0). In between around a; = one observes the crossover distribution, 
which is exactly what have been computed in this paper. 

Concretely, from (3.3), —h^^x,D corresponds to the free energy //3,-y,_D(a;, t) of the random 
directed polymer and hence (2.13) is translated as 



where F^^{x;X) is given in (2.17) and a and jt arc written in terms of 7 and B as a — 
{P^jDy/^ and jt = {P^jDH/2y/^ respectively. Based on (7.4), one can obtain various 




(7.4) 
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information about the the free energy of the directed polymer. For instance one sees that 
as ^ — >■ oo the macroscopic free energy scales as 0(/3^), the fluctuations of the bulk free 
energy is 0(/3^/^). To get a nontrivial distribution in the low temperature limit j3 — )■ oo, 
one has to scale the space direction as X = 0{/3^^^). The free energy distribution in 
this hmit is given by (2.22) since /3 — )■ oo implies 74 — >■ oo. This describes the crossover 
between boundary and bulk free energy fluctuation, i.e., between Gaussian and the GUE 
Tracy- Widom distribution. 



8 Conclusion 

In this paper we have considered the KPZ equation (1.1) with the half Brownian motion 
initial condition (1.4). Using the Bethe ansatz results of the one-dimensional attractive 5- 
Bose gas, we have obtained a Predholm determinant expression for the generating function 
of exponential moments of the height (Theorem 1). Thanks to this result and Proposition 
2, the compact representation of the probability distribution of the height was obtained 
(Theorem 3). We have also found an expression for the multi-point generating function 
employing an approximation proposed in [33, 34]. 

These results are expressed in terms of the deformed Airy functions (2.11) and (2.12). 
If we change these functions to the ordinary Airy function, we recover the results for the 
narrow wedge initial condition obtained in [14, 15, 16, 17]. The deformed Airy functions 
have some nice properties, which are discussed in the Appendix. They satisfy the biorthog- 
onality relation. Their time evolution by the Airy Hamiltonian are again given by the same 
functions with a parameter modified. We remark that similar relations of the multiple Her- 
mite polynomials played an important role in the study of the PNG model with external 
source [42, 43]. In the long time limit, the kernel of the Fredholm determinant becomes 
the rank one perturbation of the Airy kernel. It is not clear whether the higher rank 
perturbations studied in [47] have corresponding finite time generalizations. 



A Properties of the deformed Airy functions 

In this appendix, we pick up a few properties of the gamma-deformed Airy functions (2.11) 
and (2.12), which are necessary for the discussions in Sec. 6. Let t <0 and H be the Airy 
Hamiltonian H = —d^/dx^ + x. We have the following relations. 

(i) The deformed Airy function representation of the propagator 

/oo 
dze~*^ Air {x + z,h,c-ht) Ai^ {y + z,b,c). (A. 1) 

-00 

Note that when we set i = in (A.l), we get the biorthogonahty relation of the deformed 
Airy function. 



J — t 



dwAirix + w, b, c)Af{y + w, b, c) = 5{x - y). (A.2) 
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Using (A.l), we easily get the following relations, 
(ii) "Time evolution" by the Airy Hamiltonian 

e*^Ai'^(a: + w,b,c)= e-*^Ai^(a; + w,b,c-bt), (A. 3) 

e*-^Air(a; + w, b, c) = e"*'"Air(x + w,b,c + bt). (A.4) 

When we set b = l/'~ft and c = —X/'-ft and take the limit '-ft — ^ c>o, these rela- 
tions (A.l), (A. 3) and (A.4) become those for Ai~{x,y) (6.24) and Ai^{x,y) (6.25), 

dze-*"Ai+ {x + z, X) Ai" {y + z,X), (A.5) 

■oo 

e^^ki-{x + w,X)^e-^'"ki-{x + w,X -t), (A.6) 
e*'^Ai+(x + X) = e-*^Ai+(x + w,X + t). (A.7) 

Proof 

(i) The factor {x\e^^\y) can be represented in terms of the Airy function, 

/oo 
rfwe-*"'Ai(w + x)Ai(w + y). (A.8) 
-oo 

Thus it is enough to show that 

/oo /»oo 
dwe~^'"kiY{x + w,b,c-bt)A.f{y + w,b,c)= / dwe'^"" ki{w + x)ki{w + y). 
■oo J —oo 

(A.9) 

In the definition (2.11) of Ai^(a, 6, c), we find that in the case c > 0, the contour integral 
Fjc/b can be replaced by R while in the case c < 0, we can divide the contour integral into 
the integral on M and the contributions from the poles at (ic + m)/b, m — 0,1,2 ■ ■ ■ of 
r{ibz + c). Thus (2.11) can be represented as 



^ /!°^ c^^e-"+^Vr {ibz + c), c> 0, 

^ dze^-^^T i^bz + c) + E^^^.o ^^e"^-^, c < 0. 



(A. 10) 



Here n is the maximum integer satisfying c + n < 0. 

Thus for the case c > 0, using (2.12) and the first relation in (A. 10), we find the left 
hand side of (A.9) is represented as 



/oo /"OO /"OO 

-oo J —oo J —oo 



27r V{-ibzi + c - bt) 



27r 



J-oo 7]R-rf J-oo r(-z6^i + c - 6t) 
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Here we deformed the contour of zi to R — it. This deformation is vahd when t <0. Note 
that in this equation above, the part Zi + Z2 + it becomes real. Thus one has 



ir.3_u.3^ T{ibZ2 + c) 



/poo 
dzi / dz2S{zi + Z2 + it)e 
_:-it J -00 T{-ihzi + c-bt) 

r roo 

= dzi dz25{zi + Z2 + it)e^{^i^+^2j/)+i(z?+.|) 

jR-it J -00 

-| /"OO pOG poo 

= _ / dw dzi (^^2e^^(^^+^^+^*)+^(^i^+^^^)+t(^?+^l). (A. 12) 

SttJ-oo J -00 J -00 

In the second equahty, the factor of the Gamma functions is ehminated because of the 
delta function. We easily find that the last expression of this equation is equal to the right 
hand side of (A. 9) from the integral representation of the Airy function, 

Ai(a;) = / dze'^''^'-. (A.13) 
J —00 

For c < 0, on the other hand, the terms coming from the pole contributions in (A. 10) 
are added. However we find that these terms do not affect the result since 



«/ — ( 



dwe-'"^ Air{x + w, b, c)e-^(^+"') 



-j^ noo t-oo ^i(^zi+i^^^+iijw-^^^y+ixzi+i^ 

— — / dw dzi- 



27ry_^ y_oo T{-ibzi + c-bt) 

dz,6 ( zi + i^-r^ + It ——^ — = = 0. (A.14) 

'--±K+t) \ b J r{-ibzi + c - bt) r(-iV) ^ ^ 

(ii) Here we only consider (A. 3) and omit the proof of (A.4) since it can be shown by using 
the same strategy as (A. 3). The left hand side is given by 

/oo 
dy{x\e'''\y)Af{y + w,b,c). (A.15) 
-00 

Substituting (A.l) into this equation, we readily obtain the right hand side of (A. 3), 

/oo poo 
dz / dye'^"" Ai^{x + z,b,c- bt)Air{y + z, 6, c)Ai^(y + w, 6, c) 
-00 J —00 



-00 

= e-*"'Ai'^(x + w,6,c-6i). (A.16) 
Here in the second equality, we used the biorthogonality relation (A. 2). 
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